在前面《练习R：正交试验极差分析、方差分析与统计图形》一文中，用R处理了正交试验的数据，画了一个交互作用图，其实效果并不是最佳的，这个图它是矩阵的形式，看上去有两个图但是信息重叠了。

所以我需要针对方差分析交互作用显著后，直接绘制出有交互作用的两因素交互图。

这里用到的是R自带的interaction.plot()函数，语法如下（重要参数）：

interaction.plot(x.factor, trace.factor, response, type = c("l", "p", "b", "o", "c"), col = 1, pch = c(1:9, 0, letters), main = " ", xlab = " ", ylab = " ", ...)

x.factor：自变量，即因素，但记住这个变量将形成x轴； trace.factor：两因素中的另一个因素； response：因变量； type：线（l）、点（p）、线和点（b）等； col：添加颜色； pch：绘制点符号类型； main：标题； xlab：x轴标题； ylab：y轴标题。

注："p" for points, "l" for lines, "b" for both points and lines, "c" for empty points joined by lines, "o" for overplotted points and lines.

继续用design2数据，其中因素tem与oxy交互作用显著，接下来绘制交互作用图。

with(design2,interaction.plot(tem,oxy,results,col=c("black","blue"),pch=c(16,18),type="b",lwd=3))

看结果：

两条线交叉，也辅助说明前面交互作用是显著的。

本文完 文/图=数据小兵